Superconductivity in Ferromagnetic RuSr2GdCu208 



O 

o 
o 

(N 

-t— > 
o 

O 



Jian-Xin Zhu/^) C. S. Ting/i) and C. W. Chu^^'^) 
Texas Center for Superconductivity and Department of Physics, University of Houston, Houston, Texas 77204 
''^^ Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720 

The phase diagram of temperature versus exchange field is obtained within a BCS model for d-wave 
superconductivity in Cu02 layers which is coupled to ferromagnetic Ru02 layers in RuSr2GdCu2 08. 
It is found that the Fulde-Ferrell-Larkin-Ovchinnikov state is very sensitive to the band filling factor. 
For strong exchange field, we point out that superconductivity could only exist in the interfaces 
between ferromagnetic domains. The magnetization curve is calculated and its comparison with 
experiment is discussed. We also propose the measuring of tunneling conductance near a single 
unitary impurity to detect the strength of the exchange interaction. 

PACS numbers: 74.25.Ha, 75.30. Et, 75.25.Jb 
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The problem of coexistence of superconductivity (SC) 
and ferromagnetism (FM) has attracted keen interest 
since the original works of Ginzburg and Matthias et 
al. iQ. It was shown that singlet SC and FM are mutually 
exclusive and the SC can also be strongly suppressed by 
magnetic impurities. The competition between SC and 
FM were observed in the ternary compounds, HoMogSs, 
HoMogSes, and ErRh4B4 But true microscopic coex- 
istence was found only over a narrow temperature region 
when FM sets in and modifies itself to a spiral or domain- 
like structure. The recent discovery of SC (Tc = 16-47 
K) in the ferromagnetic [Tm = 132 K) ruthenate-cuprate 
layered compound RuSraGdCuaOg (Ru-1212) |-|l§ re- 
newed the interest in the issue of how SC and FM nego- 
tiate to coexist. Recent band structure calculation per- 
formed by Pickett et al. showed that the exchange 
splitting in the Cu02 layer is small (Aexc = 25 meV) 
compared to ~ 1 eV in the Ru02 layer but is larger 
enough that the superconducting state in the Cu02 layer 
may be of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
type p^Jl^ or finite momentum pairing state. Whether 
there actually exists such a new pairing state even in early 
found magnetic superconductors remains to be controver- 
sial ijl^ . Indeed, until now no superconductor has been 
discovered to be a finite-momentum pairing state. Also 
notice that, in earlier studies, the FFLO state was dis- 
cussed by assuming a constant density of states (DOS), 
that is, the exchange field does not introduce structures 
in the electronic properties on the energy scales relevant 
to SC. However, this becomes not true when the DOS has 
singular structure. Thus a further study on the existence 
of this state by including the energy dependence of DOS 
near the Fermi surface should be interesting. Moreover, 
it was supposed that the SC mainly occurs in the Cu02 
layers in Ru-1212, the identification of the pairing sym- 
metry in this material is also of fundmental importance 
in view of the well-established d-wave pairing symmetry 
in high- Tc cuprate superconductors, which have similar 
crystal structures. 

The purpose of this paper is three fold: (i) By as- 
suming a d-wave pairing symmetry in a two-dimensional 



(2D) lattice model, we present a detailed study of the 
temperature-exchange field phase diagram to invetigate 
the sensitivity of the FFLO state by varying the position 
of the Fermi energy within the tight-binding band; (ii) 
Arguing that the mixed state is intrinsic, we calculate 
the magnetization as a function of an applied magnetic 
field and compare it with experiment; (iii) We propose 
to measure the existence of the zero-energy peak (ZEP) 
and its splitting in the differential tunneling conductance 
near a single unitary impurity in the Cu02 layer as a test 
of the d-wave pairing symmetry as well as the strength 
of the exchange field. 

Our model system is defined on a 2D lattice with pair- 
ing interaction taking place between two electrons on the 
nearest-neighbor sites, which in the mean field approxi- 
mation leads to the Bogoliubov-de-Gennes equations 
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Ui6ij, and the self-consistent condition 



with the single particle Hamiltonian -ffy.o- 
fiSij — ahexcSi 

= ^E„,.K>"; tanh(£;„/2fcBT). Here 
(wicr, Viff) are the Bogoliubov quasiparticle amplitudes on 
the z-th site; 7 — zkx, ±y represents the relative posi- 
tion of sites nearest-neighboring to the i-th site; t is the 
effective hopping integral between two nearest-neighbor 
sites within the Cu02 plane; fj, is the chemical poten- 
tial; hexc — J{{S^) + {Sz)) is the exchange field coming 
from the ordered spins in the two nearest-neighboring 
ferromagnetic RUO2 layers; Ui if any accounts for the 
scattering from the impurities; V is the strength of the 
nearest-neighbor pairing interaction. Notice that the in- 
ternal magnetic field on Cu02 layers due to the magnetic 
moment on Ru02 layers is only several hundred gauss 
and the exchange interaction should be dominant in sup- 
pressing SC. We therefore defer the effect of the internal 
field on the orbital motion of paired electrons to the study 
of the magnetization. 

Temperature-exchange-field phase diagram. — For a 
pure system in the absence of external magnetic field, 
the Bogoliubov quasiparticle amplitude can be generally 
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written as 



yields the bond order parameter 
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Here N — x Ny is the number of two-dimensional 
lattice sites. Eq. (p|) shows that the order parameter is 
not a constant in space manifesting the collective motion 
of paired electrons each with momentum q. Using the 



definition = (A^ + A^"^ - A^"^ - A^"0/4, wc find 
the equation determing the d-wave energy gap 
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^k,q,- -t- -E'k.q ' '^i*^ — —2t{cOSk,j, + COS ky) — fl, 

^k,q,± = (^k+q±Ck-q)/2, and Ak = 2Ad(cosfca;-cosfcj^). 
Correspondingly, the free energy per lattice site is given 
by 
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where f^^ — /(^-k'q) Fermi distribution function. 

To determine the phase boundary between the nor- 
mal pairing (q — 0) state and the normal state (spin 
polarized), one should compare the free energies of the 
superconducting state and normal state, using Eq. (Q). 
In the presence of a fairly strong exchange interaction, 
the system might also go into the FFLO state in which 
all the Cooper pairs have a single non-vanishing (q ^ 0) 
center of mass momentum. This transition between the 
FFLO state and the normal state would be of the second 
order. To find the transition curve for the FFLO state 
and the normal state, we solve Eq. (|) with A^ = to 
find the maximum value of Ziexc by scanning through a 
whole set values of q at the same temperature T. By 
repeating the same calculation at a different value of T, 
we then obtain the phase curve /lexc as a function of T. 

To see the sensitivity of the FFLO state to the Fermi 
energy position in a 2D tight-binding band, we fix the 
pairing interaction as y = 2t and consider three typi- 
cal values of the chemical potential, fj, = —t, — 0.5t, and 



— 0.14t, corresponding to the band-filling factor v k, 0.65, 
0.82, and 0.95 (i/ = 1 is a half-filled band). For the above 
sets of parameters, the maximum energy gap at zero tem- 
perature and /lexc = is found to be Ap = 4A(i ss 0.65t, 
0.88t, and 0.96i, and correspondingly, the transition tem- 
perature is Tco ~ 0.26i, 0.37i, and 0.40t, respectively. 
Figure |l] plots the temperature-exchange-ficld phase di- 
agram. Our calculation shows that the momentum q 
corresponding to the maximum /lexc is along the (10) 
and its equivalent direction in the whole temperature re- 
gion because the energy gap reaches the maximum value 
along these directions so that the system can be more 
robust against the depairing effect from both the finite 
momentum and the exchange field. As shown in Fig. 0(a), 
when /i — —t, the transition curve (solid line) is be- 
tween the superconducting state (with q = and q 7^ 0) 
and the normal state, which shows that at low temper- 
atures, when the exchange field is increased the system 
initially in the normal pairing state will enter the FFLO 
state (q ^ 0) through the first-order transition and then 
pass into the normal state by a second-order transition. 
Thus the FFLO state is a stable state at high exchange 
fields. The transition curve between the normal pair- 
ing state and the FFLO state is represented by dashed 
line. When the temperature is increased, two curves be- 
comes closer and at T = 0.58Tco (tri-critical point) they 
coincides with each other, where the transition begins 
to be of second order and the FFLO state merges nat- 
urally to the q = normal pairing state. This result 
is similar to that obtained within the continuum model 
using a constant DOS |T2|jr^ ]. However, wc find that the 
phase space for the existence of the FFLO state shrinks 
ai ji — —Q.ht as the band-filling factor shifts toward the 
half filling (Fig. |l](b)). In particular, near half filling 
when ji = — 0.14i, the FFLO state will be unstable and 
appear as a supercooling state. Therefore, as shown in 
Fig. 0(c) , only the transition between the normal pairing 
state and the normal state is physically acceptable. The 
novel feature comes from the influence of the exchange 
field on the DOS near the Fermi surface. For a 2D tight- 
binding band, the DOS at the energy /i = has a singu- 
lar point and it decays logarithmically as the energy goes 
away from zero. When /i = — i, which has been far away 
from the zero energy point, the DOS is flat for the en- 
ergy near the Fermi surface. In this case, the splitting of 
the normal electron band by the exchange field has little 
effect on the change of the DOS and a constant DOS can 
be taken which corresponds to the approximation made 
in the continuum theory (l^Jl^ . But as the chemical po- 
tential is close to the zero energy point, a little splitting 
of the energy band will cause a strong variation of the 
DOS near /i = 0, which makes the q ^ pairing state 
unfavorable. 

For the Ru-1212, the band structure calculation esti- 
mated the Fermi wavevector as tt/o. The filling factor 
should be near the half-filling (/x = 0). Then the exis- 
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tence of the FFLO state in the above system becomes un- 
hkely according to Fig. |^(c). In the following discussion, 
we focus on the transition as between the normal pairing 
state and the normal state. Consider the case of = 0.95 
(Fig. |l|(c)). If the zero- field transition temperature is as- 
sumed to be Tco — 90 K, Aq is about 18.7 meV. By fitting 
to the transition temperature ~ 36 K as measured ex- 
perimentally 1^ , the exchange field should be as small as 
8.8 meV. Actually, Tcq may be as low as 60 K when Ru is 
replaced by other atoms (e.g., Cu) ||l^. Then the allow- 
able hexc is only 5.1 meV to have ^ 36 K. Notice that 
the exchange field hexc from the band structure calcula- 
tion is as large as 12.5 meV |11 and the experimental 
filling is very close to = 1 (or = 0), we therefore con- 
clude that, if hexc is indeed so large, SC can only exist 
in the Cu02 layers between the ferromagnetic domains 
where the exchange interactions hexc = J{{Sz) + (S"^)) 
are small. It is important to emphasize that we have not 
ruled out the possibility of the FFLO state in the real sys- 
tems because the estimated value of exchange interaction 
could in fact vary over a limited range depending on the 
method of obtaining it. In case that the real exchange 
interaction is somewhat smaller than that estimated in 
Ref. [0 , the bulk FFLO state would become a reality. 

Magnetization. — The bulk Meissner effect has not 
been observed in the superconducting state of Ru-1212 ||] 
until very recently ||^. Furthermore, even the absence 
of superconductivity has been reported in well charac- 
terized Ru-1212 samples |l^. This discrepancy indi- 
cates the delicate balance between the superconduct- 
ing and ferromagnetic interactions and the experimen- 
tal result appears to depend critically on the sample 
condition in a yet-to-be determined fashion. It is re- 
ported that the magnetic moment in each Ru atom is 
1 Bohr magneton [Q. With structure parameter values, 
a = 3.8 A, and c = 11.4 A, the internal magnetic field 
Hint = 4:TrMsp = l/is/a^c {Mgp is the spontaneous mag- 
netization) is estimated to be 707 G in the Cu02 layers 
which is larger than the first critical field H^^'' ~ 100 G of 
a non-Ru layered cuprate superconductor with the com- 
parable transition temperature, i.e., Hj,^' — AnMsp < 0. 
The measured Hd is therefore zero and no bulk Meiss- 
ner effect can be observed. Instead the superconductiv- 
ity occurring in the Cu02 layers has been driven into the 
mixed state. The overall magnetization in the system 
consists of two parts, one from the spontaneous mag- 
netization Msp of the ferromagnetic Ru02 layers, the 
other Mob from the diamagnetic orbital contribution of 
the superconducting Cu02 layers in the mixed state, i.e., 
M = Msp + Mob- When an external magnetic field Hext 
is applied, the effective magnetic field is H — H^xt+Hint- 
As an approximation, we work with the London equation 
for a square vortex lattice to find the magnetic induc- 
tion B K. H - Hciln{Hc2/B)/lnK 0, where k is the 
Ginzburg-Landau parameter and Hc2 — 2k^Hci is the 



upper critical field for the Cu02 subsystem. In Fig. the 
total magnetization is plotted as a function of Hgxt- As 
it is shown, the magnetization increases monotonically 
with the external magnetic field. The observed mono- 
tonic behavior of M by the experiment Q is shown in 
the inset of Fig. ||. Because our study is based upon the 
SC in a single ferromagnetic domain, the internal mag- 
netic field must be greater than the intrinsic first critical 
field Hci of the Cu02 SC. Nevertheless, when Hext — 0, 
our calculation gives an appreciable spontaneous magne- 
tization which contradicts the experimentally measured 
zero magnetization This difference shows the exis- 
tence of ferromagnetic domains in Ru-1212, where the 
average magnetization vanishes in the absence of an ex- 
ternal magnetic field and the SC occurs in the interfa- 
cial regions between some of the ferromagnetic domains. 
Note that the magnetization has also been recently stud- 
ied |l^ to analyze the superconducting properties of 
Ri.5Ceo.5RuSr2Cu20io, where the Meissner effect was 
absent at temperature region Td < T < Tc (Tc ^ 30K 
and Td ^ 20K) but present at T < Td- To interpret this 
phenomenon, it was proposed that the ferromagnetism 
appears with a domain structure but the superconduc- 
tivity is a bulk phase, which is different from our expla- 
nation for Ru-1212 systems, where neither Meissner effect 
(down to 0.5 Gauss) nor detectable condensation energy 
was observed at temperature down to 2K Very re- 
cently, the detection of Meissner state at a field below 30 
Oe was reported by Bernhard et al. In the reported 
sample with a Meissner effect, the internal field was esti- 
mated to be only about 50-70 Oe (the lower critical field 
Hci of the nonmagnetic superconductor was estimated to 
be of the order 80-120 Oe), in contrast to the previous 
reported Hint about 200-700 Oe by the same group and 
others. We argue that the experimental result of Bern- 
hard et al. might make sense if the superconductivity 
occurs in the inter-ferromagnetic domain region where 
Hint is much reduced and the ratio R between supercon- 
ducting volume/sample volume is not too small. This 
would not be inconsistent with the conclusion reached 
for Ref. m, where R could be rather small so that the 
Meissner effect was not detected. 

Quasiparticle resonant state near a single unitary im- 
purity in Cu02 layers. — Since the ferromagnetic ex- 
change field has pre-existed in the above system, it can 
affect the quasiparticle resonant states near a single im- 
purity in the case of d-wave pairing symmetry. To ad- 
dress both issues, we solve Eq. (^ self-consistently using 
the exact diagonalization method ||l^ and calculate the 
local density of states (LDOS) Pi = ~ Yl,n cr[l"jo-P/'(£^ — 
En) + \vis\'^ f'{E + En)]. The calculation was made on 
6x6 supercells each with size 35a x 35a by assuming 
a paramagnetic pairing state. The parameters are as 
follows: The single-site impurity strength Uq = lOOi, 
V = 0.95, and T = 0.02t. From Eq. (0), we can see 
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that the zero-field quasiparticle energy E^^^ is shifted to 
hexc- Therefore, the position of zero-energy 
states is now spUt to ±hexc- Figure^ plots the LDOS on 
the site nearest-neighboring to the impurity site. As is 
shown, when /lexc — 0, there occurs a single ZEP in the 
LDOS. In the presence of exchange field, the ZEP is split 
into double peaks, each corresponding to one spin compo- 
nent. The magnitude of the splitting increases with the 
exchange field. Since the local differential tunneling con- 
ductance is proportional to pi, the ZEP and its splitting 
can be detected by the scanning tunneling spectroscopy 
(STS) as a test of the d-wave symmetry as well as the 
pre-existing exchange field in the above system. Experi- 
mentally, the nonmagnetic impurity with strong scatter- 
ing potential can be realized by substitution of Zn for Cu 
in the Cu02 layers. In addition, the STS best suited for 
exploring the local electronic properties allows a direct 
examination of whether the SC in Ru-1212 appears as 
a bulk state or can only survive at the boundaries be- 
tween ferromagnetic domains. For the former, the STS 
data should reveal a superconducting gap over the whole 
sample. 

Finally, we would like to mention once again that the 
bulk SC in Ru-1212 prevails only when the exchange in- 
teraction is weak. For large exchange interaction, the SC 
could only exist at the interfacial Cu02 layers between 
ferromagnetic domains, where /lexc is small. Whether 
the SC in this compound is bulk or interfacial like ap- 
pears to depend on sample preparation and this issue 
needs further experimental studies. 

Note added: After the paper was submitted for publi- 
cation, we noticed a preprint by Shimahara and Hata pC| ] 
in which the enhancement of the possible FFLO state in 
a layered ferromagnetic compound such as Ru-1212 by 
the next-nearest-neighbor hopping was discussed. 
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FIG. 1. The T-hexc phase diagram for a d-wave supercon- 
ductor with fi = -t (a), -0.5t (b), and -0.14i (c). The 
solid line represents the physical phase boundary between the 
superconducting state and the normal state. In panels (a) 
and (b), the dashed line indicates the transition between the 
FFLO state and the normal pairing (q — 0) state. Note that 
the FFLO state exists only as the temperature is below the 
tri-critical point A. In panel (c), only the transition between 
the normal pairing state and the normal state is allowed. 



FIG. 2. Magnetization as a function of the external 
magnetic field. The parameter values: Hci = 100 G, 
Hint = 707 G, and k — 50. The inset shows the measured 
magnetization |8|. 



FIG. 3. Local density of states as a function of energy for 
hexc ~ (solid line), O.lAo (dashed line), 0.2Ao (dotted fine) 
at the site one lattice constant away from the impurity site. 
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FIG. 1 Zhu et al. 




FIG. 2 Zhu et al. 




FIG. 3 Zhu et al. 



